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The entanglement between two arbitrary subsystems of random pure states is 

rp 

studied via properties of the density matrix's partial transpose, p^2- The density of 
states of pj'2 is close to the semicircle law when both subsystems have dimensions 
which are not too small and are of the same order. A simple random matrix model 
for the partial transpose is found to capture the entanglement properties well, in- 
cluding a transition across a critical dimension. Log-negativity is used to quantify 
entanglement between subsystems and approximate analytic formulas for this are 
derived. The skewness of the eigenvalue density of is derived analytically, using 
the average of the third moment that is also shown to be related to a generaliza- 
tion of the Kempe invariant. Extreme value statistics, especially the Tracy- Widom 
distribution, is found to be useful in calculating the fraction of entangled states at 
critical dimensions. These results are tested in a quantum dynamical system of three 
coupled standard maps, where one finds that if the parameters represent a strongly 
chaotic system, the results are close to those of random states, although there are 
some systematic deviations at critical dimensions. 
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I. INTRODUCTION 

Quantum entanglement is a central property of quantum mechanics that is absent in 
classical physics. Studied since Schrodinger and the famous paper of Einstein, Podolsky 



and Rosen (EPR) jl|, correlations due to entanglement seem to imply nonlocality. The 
work of Bell and others led to inequalities that quantified the extent to which classical 
correlations can be surpassed. These inequalities were experimentally verified by Aspect 



et.al [Sj. However, entanglement has been extensively studied more recently as it is a critical 
resource for quantum computation ^, quantum teleportation js], dense coding jsl, and 
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maximally entangled state of two qubits. 



various other_^quantum information tasks [71, Isj, and to explain the magnetic properties of 
some solids [9.]. A well known example of an entangled state is the spin singlet which is a 



It is known that in a generic or random pure 



state any of its subsystems is nearly maxi- 



mally entangled with the complementary system |l0l4l3l|. where the measure of entanglement 
is the von Neumann entropy of the subsystem. Here "random" means that the state is sam- 
pled uniformly from the unique Haar measure that is invariant under unitary transforma- 
tions. Random states occur in many contexts. For exanaple, they are found as eigenstates 



of quantum maps whose classical limit is fully chaotic [IJ]. For the eigenstates of quan- 
tum systems with classically chaotic, continuous Hamiltonian analogs, one must account for 
an effective dimensionality, i.e. energy window such as the Thouless energy, in addition. 



With that proviso, disordered or chaotic ballistic mesoscopic systems [15| exhibit random- 
ness in their single particle eigenstates, and in strongly interacting systems such as medium 
to heavy nuclei with many valence nucleons, there is randomness in the full many-body 
eigenstates [3]. 

The interest in this paper is to study the entanglement between two blocks comprising say 
Li and L2 qubits in a random pure state {ip) of L qubits (Li + L2 < L) (see Fig. ([T])). While 
one can work with any dimensional Hilbert space, this paper puts the results mostly in terms 
of collections of "qubits" or spin-1/2 particles, the generalizations being straightforward. 
The reduced density matrix, pi2, of Li + L2 qubits is obtained by tracing out the remaining 
L — Li — L2 qubits: 

Pl2=tTL-L,~L,m{i;\). (1) 

The state pu is in general a mixed one, i.e. tr(pi2)^ < tr(pi2) = 1, and the equality holds 
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only if the qubits in the blocks 1 and 2 are unentangled from the rest. A mixed state of a 
bipartite system is separable if and only if it can be written as 

Pi^ ® pf\ with Pi > and = 1, (2) 

i i 

(1 2) 

where p] ' are density matrices of subsystems 1 and 2. Otherwise the state is non-separable, 
or entangled. Given a general state it is a challenging task to verify if it is separable or not. 

One simple (partial) test for entanglement is Peres's partial transpose (PT) criterion jirj l. 
The matrix transpose map T : p — )■ is trace preserving and positive, since for every p > 0, 

> 0. However its extension / T to a bipartite system (where I is an identity matrix 
that acts on the first subsystem and T acts on the second subsystem) does not preserve 
positivity. Hence transposition is a positive but not a completely positive map, and can be 
used to reveal entanglement. The map / T is called a partial transposition (PT) since it 
effects transposition only on the second subsystem keeping the first subsystem unaltered. 
The test is partial as it leads to necessary but not sufficient conditions for entanglement. 

Given a bipartite system 1 and 2 having an orthonormal basis {|^)|q;)} and density matrix 
Pi2, the PT with respect to the second subsystem, denoted as pfa, is given by the matrix 
elements: 

(pS)i«;i/3 = (Pl2)i/3;ja ; ipl2)ia;j(3 = | (« |pi2 1 j) | /?) • (3) 

Peres's partial transpose (PT) criterion states that if p§ is negative then the state pu is 
entangled. In this case pi2 is said to be a NPT (negative partial transpose) state, otherwise 
it is a PPT (positive partial transpose) state and is guaranteed to be separable only for 2 x 2 



and 2x3 systems 



concurrence 
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Lsjl8|. 



Entanglement between 2 qubits in a mixed state is also given by the 
20| which takes values from to 1, where corresponds to an unentangled 
or a product state and 1 corresponds to a maximally entangled state. For more than two 
qubits or higher dimensional quantum spins in a mixed state, negativity and log-negativity 



21 



221 ] are used as measures of entanglement. 
In the following, eigenvalues of a density matrix without PT are denoted by Aj and those 
of pfg by Pi. Negativity of the state pi2 is defined as the sum of the moduli of the negative 
eigenvalues of p^2 5 which is clearly zero for PPT states. Due to the fact that the trace is 
preserved under partial transpose the negativity is also 

^fip^2) = (4) 
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Log-negativity is defined as 

Ei;v = log(||pr^||i) =log(^|/i,,|). (5) 

i 

If tlie log-negativity is greater tlian tlien tlie density matrix is entangled. Otlierwise tlie 
state pi2 is separable or it could also be bound entangled |23|. Bound entangled states are 
entangled but they can not be distilled by means of local operations and classical commu- 
nication to form a maximally entangled state. The distribution of the eigenvalues of 
of central concern in this paper. Note that the trace of the first and second powers of the 
density matrix remains unaltered under the PT operation. The first power to show deviation 
between the two sets of eigenvalues is the trace of the third power. That is: 

tr (pi2) = tr [pID = 1, tT{p,^) = tT[{pJlf], tripun + ^A^pllT] m>3. 

The average of tr[(p^2)^] is explicitly evaluated further below, for both real and complex 
states, where the average is over all the pure states (see Eq. ([1])) sampled uniformly. 
Interestingly this quantity is a generalization of one of local unitary invariants studied for 
three qubits [24], and therefore is of broader interest. 

The distribution of the eigenvalues, i.e. the density of pfg (that of pj), is of evident 
importance in a calculation of the entanglement between subsystems 1 and 2. In this paper 
a simple random matrix model is proposed for the partial transpose, based on the known 
average of the second moment. This model quite accurately predicts a transition from 
dominantly NPT states to dominantly PPT states as the size of the subspaces Li and L2 
are varied. The transition region is an interesting one wherein the extreme eigenvalues 
of random matrices determine the nature of the entanglement. Use is made of the well- 
known Tracy- Widom distribution to estimate the fraction of NPT states in the transition 
to predominantly PPT ones. The limitations of the simple model are also pointed out, 
especially when the skewness of the densities are important and Li and L2 differ significantly. 

Finally in this paper a dynamical model of three coupled standard maps or rotors is stud- 
ied, restricting attention to the case when they are classically fully chaotic. The eigenstates 
of such a system are taken to be the pure states in Eq. ([T]) and the entanglement between 
rotors is studied via the log-negativity measure. While good agreement is found away from 
the transition region, there are interesting deviations in this critical zone. While all standard 
diagnostics, such as the distribution of the nearest neighbor spacings of the eigenangles, the 
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y. Recently Datta in 26 1 



number variance, the distribution of the eigenvector components, agree with random matrix 
theory (RMT) to a large extent, deviations are seen with respect to the fraction of NPT 
states. Stated simply the dynamical system has systematically more entanglement than 
predicted by random matrix theory. These tests are perhaps some of the more stringent 
ones of the Bohigas-Giannoni-Schmit conjecture {25] that quantized chaotic systems have 
spectra whose statistical properties are modeled by those of random matrices. These tests 
are stringent as they rely on outliers or extreme eigenvalues. In the large system dimension 
limit (small effective h) there does, however, seem to be a tendency to approach the RMT 
results. 

Works related to the results in this paper have appeared previous 
has studied entanglement of random pure states using negativity 
which adds to the full system (Li = L2 = L/2) and found that the average negativity 
is a constant (0.72037) multiple of the maximum possible (= (2^/2 _ 1^2). This IS a 
reflection of the large entanglement present in random bipartite pure states, for example 
as measured by the von Neumann entropy of the subsystems A calculation presented 
further below, based on previously derived results in 27| , gives an explicit expression for the 
average negativity that is also slightly different. In this case the eigenvalues /ij of the partial 
transpose are simply related to the eigenvalues of the reduced density matrix and therefore, 
imphcitly, the density of the eigenvalues /ij had been derived even earlier |28|._ 

If the two subsystems do not make up the full system, Kendon et. al. in 



17( 1 for equal bipartition 
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30| found 



numerically that in a typical random pure state the subsystem consisting of Li and L2 qubits 
is NPT if Li + L2 > -L/2. Analytically they showed that the lower bound on Li + L2 for 
P12 to be NPT is L/3. It is shown in this paper that using the simple random matrix model 
for the partial transpose leads to the bound on Li + L2 for pi2 to be NPT is in fact L/2. 
In 31] Carteret has given a quantum circuit which can determine the spectrum of p§ by 
computing tr(p^|)' for all /'s up to the dimension of p§. Then from Peres's partial transpose 
criterion one can determine whether pi2 is NPT or PPT. 

The random mixed states studied in this paper are those that arise from a partial trace 
of random pure states selected according to the Haar measure. Properties of random mixed 
states generated according to the measure induced by the Bures metric [32] , have been stud- 



ied earlier using the von Neumann entropy and purity in 32, l33| . Multipartite entanglement 



for localized states 



3J] (35|, and multifractal states (using the von Neumann entropy) have 
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also been studied. 



. Mathematical work connected 
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:o the spectrum of the partial trans- 
which is of a complementary 



pose has appeared very recently in the literature 
nature, but with some overlap, after much of the present work was done. 

The structure of the paper is as follows. In section [Tll some known and relevant results on 
the reduced density matrix are first summarized. Rest of this section is a detailed treatment 
of the effect of PT on the reduced density matrix, in particular a random matrix model is 
seen to give rise to the observed Wigner semicircle density of states on PT, and predicts 
the transition from a predominantly NPT to primarily PPT phase. Further in this same 
section we calculate the average of the the trace of the third power of the density matrix 
after PT, and show how it is related to an invariant, the Kempe invariant, that has been 
studied earlier in the literature. In section IIIII these results are used to find the average 
log-negativity between two subsystems of the tripartite state. In subsection IIII Bl results on 
extreme value statistics of minima of reduced density matrices after PT are presented, and 
it is seen how the Tracy- Widom distribution gives rise to the fraction of NPT/PPT states at 
critical dimensions. In section IIVI we compare our results of random states with eigenstates 
of three coupled quantum standard maps, and find good agreement. 



II. STATISTICAL PROPERTIES OF THE PARTIAL TRANSPOSE 



A. On the reduced density matrix of a subsystem 



If a bipartite quantum system of Hilbert space dimension N x M {N < M) is drawn 
from the ensemble of random pure states then the joint probability density function of the 



eigenvalues 
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40| of the reduced density matrix p^v of a subsystem of dimension is 



^ N N 
^ i=l ^ i=l 



|(M-Af+1)-1 



(6) 



where /3=1, 2 and 4 for real, complex and symplectic case respectively. The normalization 



constant Km,n is calculated using Selberg's integral |40|. The density of the eigenvalues, for 
large N and M, is given by an appropriately scaled Marcenko-Pastur (MP) function 271 141|. 

iVgv/(A+-A)(A-A_) 

(7) 



/(A) 



2tt 



N 



A± = TT 1 + TT± 



Q VQ 
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12 3 L-2 L-1 L 

FIG. 1. System of L qubits. Li + L2 < L 



where A G [A_,A+], Q = M/N and Nf{\)dX is the number of eigenvalues in the range A 
to A + dX. For Q = 1 there is a divergence at the origin. For Q > 1 the eigenvalues are 
bounded away from zero. 

The purity of the subsystem density matrix pat is always larger than and less than 1. 
The minimum value being attained when Pat is maximally mixed, and the maximum when 
the subsystems are unentangled. The average purity of the subsystem Pa? is given by 10 1 

The subsystem entropy is a good measure of bipartite pure state entanglement and remark- 



ab 
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y there is an exact formula for its average evaluated over the probability density in Eq. ([6]) 



42 



43|. 

^ 1 N-1 N 
{-tripN log Pn)) = E - - ^ log(^) - ^ ^ « ^ ^ ^- 

m=M+l 

In terms of an interpretation, there is practically very little information about the full pure 
state in a subsystem, to be more precise there is less than one-half unit of information on 
average in the smaller subsystem of a total system in a random pure state. The maximum 
entanglement being log(iV), there is near maximal entanglement between any two subsystems 
of a random state. 



B. Effect of PT On Reduced Density Matrix 

Reverting back to the notation of pi2 as the reduced density matrix, while its density 
of eigenvalues is the scaled Marcenko-Pastur distribution in Eq. ((Tj), we are interested in 
the spectrum of its partial transpose, pfg- It is numerically found that for Li = L2, the 
eigenvalue density of fits the well-known Wigner's semicircle law for any L such that 
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FIG. 2. Density of states of pi2 (left) and p^l (right) for Li = L2 = 3, and where L is the total 
number of qubits. A vertical line at the origin has been shown in the right figure to draw attention 
to the negative part of the spectrum. 

L1 + L2 ^ L; see Figs. Oscillations are found about the semicircle for very small values 

of Li, just as in the case of the canonical ensembles of RMT 4J, |45|. The figure shows results 
for Li = L2 = 3 and varying L from 12 to 16. The semicircle's are fit according to a center 
(or shift) and width that is discussed further below. The rather good agreement with the 
semicircle law for the spectrum of the partial transpose is evident. For instance, in the case 
when Li = L2 = L/A (corresponding to the case L = 12 in Fig. (EJ) the rescaled eigenvalues 
X = iV/i fit the following formula: 



P^[x) = — ^/4- (x- 1)2 where - l<x<3 andA^ = 2^'+^^. 
2ti 



(10) 



Recently Aubrun 371] has used the binary correlation method |16i] to find a shifted 
Wigner's semicircle law under PT. We however use an approximate and simple model that 
enables us to see the transition that is observed in Fig. ([2])) when the total number of qubits 
is L = 14. The said transition is from a predominantly PPT phase when Li+ L2 < L/2 — 1 
to a predominantly NPT one when L1 + L2 > L/2 — 1. The critical case is an interesting one 
that is fit for the application of extreme value statistics to find the fraction of NPT states. 
In this case the semicircle lower bound is at 0. When L is odd however one finds that there 
is no Li + L2 which is critical in this sense; instead for Li + L2 < {L — 3)/2 the states are 
predominantly PPT and if Li + L2 > (L — 3)/2 are predominantly NPT. If one is given a 
certain number of qubits Li + L2, then there is always the case when the total number of 
qubits is L = 2(Li + L2) + 2 which is critical. In this work most of the calculations are for 
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L even, and there is a critical subspace dimension Li + L2=L/2 — 1. When Li + L2 = so 
that the "subsystem" 1 + 2 is, in fact, the whole system and is in a pure state, much can be 
said about the spectrum of the partial transpose. This case, discussed later in this paper, 
has a density of states that is not the Wigner semicircle. However, a semicircle is obtained 
even from small deviations of Li + L2 away from L. 



1. Degree of partial derangement in the partial transpose 

The PT operation partially rearranges the positive matrix pi2 through selective exchange 
of matrix elements. One may expect that the extent of such a rearrangement will be con- 
nected with a deviation from the Marcenko-Pastur distribution and approach toward the 
semicircle law. In other words, the number of elements exchanged by the PT operation 
results in a loss of the particular correlation among matrix elements necessary to make the 
original matrix positive. However, the rearrangements do preserve the Hermitian nature of 
the matrices. Additionally, for a density matrix of M qubits, the eigenvalues of the matrix 
obtained after PT on k qubits are the same as after doing the partial transpose on the 
complementary M — k qubits. Thus, the range < A; < M/2 is the full range of interest. 

First divide the whole matrix of dimension 2^^ x 2*^ into matrices of dimension 2^ x 2^; 
the number of such matrices being 2'^^'^^^^\ PT on k qubits is a full transpose on these 
2k ^ 2'^ matrices. Therefore the number of elements getting exchanged after PT is 

n ^ ^2M _ c.2{M-k)c.k ^ _ (-^-^^ 



This number, which is evidently the same for whether the density matrix is real or complex, 
is maximum when k = M/2, and therefore one can expect the maximum loss of correlation 
among matrix elements of pi2 and the development of the Wigner semicircle law. When k 
is smaller one still obtains qualitatively different spectra depending on the density matrix. 
In Fig. ([3]) this is seen, with M = Li + L2 = 8 and k = L2. As the number of qubits k 
varies from 1 to 4. The obtained densities on PT are all very similar except for the extreme 
case of k = L2 = 1, when the skewness is more apparent. Remarkably the minimum of the 
distributions remain unchanged even as the maximum shifts slightly. The question of the 
skewness is addressed further below. 
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Wigner' s Semicircle 




FIG. 3. Density of states after partial transpose for various Li, L2 and fixed L. The skewness is 
minimum for Li = L2 = 4 and maximum when Li = 1 and L2 = 7. Except for the case when 
Li = 1 and L2 = 7 all other cases are close to Wigner's semicircle law. 

2. A model for the shifted semicircles 

A simple model for the spectral density of p^g? of a density matrix, is suggested 

by the fact that the first two moments do not change under the operation of PT. As a 
semicircular density depends on just the two moments of mean and variance, it is proposed 
to shift and scale the semicircle of the Gaussian random ensembles to match the first two 
moments of pi2 ( or equivalently p§). In particular, we assume that these random matrices 
belong to the Gaussian Unitary Ensemble (GUE). Thus consider 

B = A+^-^ (12) 

where A is a N x N GUE random matrix with the necessary matrix element variance to 
match the variance of pu and 1^ is the identity matrix. It follows that (tr(i?)) = 1 since 
(tr(yl)) = 0, where the angular brackets indicate the ensemble average. The fact that tr(i?) 
is not exactly equal to unity for each and every member of the ensemble would be expected 
to have an influence only in the case of very small dimensional cases. 

As the eigenvalues of B are all those of A shifted by it is sufficient to consider the 
spectrum of A. Under the assumption that it is from the GUE it follows that the density of 
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eigenvalues of B is (for large A^) 



where 



R = mj^MA^)) = ^Jl^{tM2)) - (14) 



Use is now made of the approximate form of the average purity in Eq. ([8]) to derive that 

(tr(A2)) = R= = 2-^/2+1, (15) 

where one has also used that = 2^'-^, N = 2^'^^^^ and Li + L2 + L3 = L. If the scaled 
variable x = fiN is used, the resultant semicircular probability density has a shift of 1 and 
a rescaled "radius" R = NR = 2^i+^2-{W2-i). Exphcitly: 



Pr(x) = -^Jr^ -{x- 1)2, l-^<x<l + ^. (16) 

This is the Wigner semicircle law that has been used in Figs. dH [3]) and illustrates how well 
this simple model works. 

Moreover this treatment gives the PPT - NPT transition as well. For if L is an even 
integer and L1 + L2 = L/2 — 1 then R = 1 and the radius of the (rescaled) semicircle is 
such that the lower limit is exactly at 0. For any L1 + L2 > L/2 — 1 the radius is larger 
than unity and there are NPT states, while in the opposite case the lower bound is such 
that there are predominantly PPT states. Thus the transition is clearly indicated in the 
model of the partial transpose as a shifted random matrix of the GUE kind. If L is odd, 
it is clear that there are no Li, L2 such that the radius is unity, but it is easy to find that 
when Li + L2 = [L — l)/2, the radius is \/2 and hence the states are predominantly NPT, 
while when Li + L2 = {L — 3)/2, the radius is 1/a/2 and hence the states are predominantly 
PPT. These are precisely the values found from numerical simulations, and again the simple 
model predicts them correctly. 

An additional interesting feature is that the model predicts that the range of the eigen- 
values is the same both before and after the PT. Namely 



N{\+ -\_) = 2R = A^/nJWs = 2^i+^2-(^/2-2) ^ (17) 

where X± are the limits of the Marcenko-Pastur distribution in Eq. ([7]). This is borne out in 
Fig. ([2]). While this is not an exact equality, it seems to be nearly true statistically. Extreme 
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deviations from this will occur when the subsystem 1 + 2 is nearly pure or pure, a case we will 
discuss later. Is there some characteristic of the density matrix pi2 that signals the PPT - 
NPT transition? Note that when Li + L2 = L/2, pu has a density of states that diverges 
at 0, see Fig. ([2]), and for Li + L2 > L/2, the density matrix is rank deficient. Whereas the 
critical case as far as this transition goes is at Li + L2 = -L/2 — 1 when the density of states 
of P12 does not diverge at zero. While the rank of the density matrix mattered in the case 



of an entanglement transition observed for definite-particle states recently [46|], it seems to 
be not exactly the case here, as there is a case when the density of states of pu is bounded 
away from zero, but its partial transpose has a significant measure of negative eigenvalues 
and is predominantly NPT. 

As is apparent from the Fig. ([3]) the semicircle is not obtained when one of the subspaces 
is of very low dimensions, although interestingly even in this case the minimum eigenvalues 
after PT remains nearly the same. We limit most of our discussions to those cases where 
the semicircle law is approximately valid. Another instance where the semicircle law is not 
valid is when the third subspace has no qubits, that is the state pi2 is itself pure. This case 
will be discussed in the next section. More work needs to do be done in elucidating the 
boundaries of the applicability of various densities after PT. 



3. The third moment, the Kempe invariant, and the skewness 

The lowest ordered moment which changes after PT is the third moment i.e. tr [(pi2)^] 7^ 
tr[(p^2)^] ^iid it is therefore interesting to calculate the exaci ensemble average (tr [(p^g)^] )) 
and compare it with that of the simple model above. In the case of complex random pure 
states of L = Li + L2 + L3 qubits, we find 

, . r,N3, _ Nf + Ar| + Nj + 3N^N2N:i _ + + 2^(^-^i-^2) +3x2^ 
^''^P^^'l ' ~ {N,N2Ns + l){mN2Ns + 2) (2^ + l)(2^ + 2) ' ^ 

where Ni = 2^\ Details of the derivation are relegated to the Appendix. In contrast, prior 
to PT, 

^''^P''^ ^ - {N,N2N, + l){N,N2N, + 2)^ ^^^^ 



so that 



(trio ? tr(o^^?) - (Nf-^ml-l) . . 

(tr(pi2) - tr(p,,) ) - (^^^^^3 + i)(^^^^^3 + 2)- ^2°) 
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Thus, on average the third moment after PT is smaller than that before. The equation in 
Eq. f|T9l) is a special case of Eq. f|T8|) . with the identification of A^2 = 1 and A^^i = N1N2, as 
the original density matrix is the same as a partial transpose over zero qubits. 

The ensemble average of the third moment after PT has a permutation symmetry as is 
clear from Eq. ( |T8|) . Quite remarkably, this is true for every realization in the ensemble, and 
is a property therefore of pure states split in a tripartite way. To be explicit, in this case 
the following can be shown to be true: 

tr(pr|)^ = tr(pg)^ = tr(prO' (21) 

Note that there is no such constraint for the density matrices pu, P23, and pis themselves. 
Numerical calculations, as well as special cases, indicate that this is a peculiar property 
of the partial transpose and only for the third moment, apart from the trivial case of the 
first moment. To our knowledge, this has been identified as one local unitary invariant 

n 

for the case of three qubits |31|, but not for general tripartite systems. In the case of three 



qubits this quantity, which has however 



invariant" and denoted as h \2A 



Deen written differently, has been called the "Kempe 



48|. 



For completeness a proof is now supplied for the identity in Eq. ( 12T1) . Let the pure 
tripartite state and its adjoint be written in a standard basis as 

\i') = Y.'^M3kl). (^| = $^?''(jA:/|, (22) 

jkl jkl 

where 1 < j < A^^i, 1 < < Ai"2, 1 < / < A^a, and if) is the complex conjugate of if). The 
following then ensues (repeated indices are summed over): 

P12 = E pg = E i^^k'iW'^'mifk'V (23) 

jkJ'k' jk,j'k' 

W T2\3 , -jj'kl I -j"k'l' -jk"l" , . 

tr(Pilj =Vjk'iip VJj'k"i'ip Wj"ki"V ■ (24) 

Similarly it follows on tracing out the second system and taking the partial transpose with 
the first that: 

. f Ti\3 , -jjki' 1 -j'k'l" -j"k"l . . 

tr(P3ij =Vj'kiV Vj"k'uV Vjk"i"V ■ (25) 

To see the equality of Eq. and Eq. (!25|) the following permutation of the dummy indices 
suffices: (j — )■ j", j" — )■ j', j' — t- j), and {k — > k', k' k", k" k). It seems somewhat 
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unusual to express the Kempe invariant in terms of the partial transpose, but this indeed 
seems to be a simple way of doing so. Since the quantity is both a local unitary invariant and 
invariant under permutation of the systems it serves as some kind of entanglement measure 
in itself. According to [47] it is a measure of bipartite entanglement, indeed it is possible 
that it is some overall measure of entanglement between any pair of the tripartite system. 
Thus the average of this quantity as found in Eq. (fT8|) is of larger interest as well. The 
average Kempe invariant of three qubits is 2/5 while for three qutrits it is 27/203. 

One may also compare with the third moment for the shifted GUE matrices in the model 
of Eq. (fT2|) . The third moment is approximately: 



1 



(tr(S^)) ^ + (26) 

\ \ ^/ AT AT AT AT2at2^ \ ' 



neglecting higher order terms. While this does have the correct leading behavior, it is not 
the same as the exact moment, and does not also possess the permutation symmetry noted 
above. The difference between this result and the exact moment is of a lower order than the 
exact moment: 

Hpa'-t..(B»))«^(^ + ^). (27) 

The skewness of the density of states is zero for the shifted GUE ensemble, but is nonzero 
for the PT of the density matrices. Skewness 7 of a distribution is the normalized third 
central moment: 

1=1 ^ ' 

where is number of elements and o is standard deviation of the sample. Using 71 = 
tr(pi2)/A^ = 1/A^, and = ii{p\^lN -JI^ ^ ^(A^iA'aA's) (recall that = N1N2) leads to 

^ ^ A2; ■ ^ ^ 

In terms of number of qubits the result is that for large Li, L2 and L 

{2-L/2^2^i-L2 _^ (complex states) 

(30) 
2-L/2^2^i-^2 + 2^2-Li + 3(2-^1 + 2-^2)] (^real states). 

The case of real states is dealt with the Appendix where further details are found. Thus it 
follows that for a given L and Li + L2, the skewness is a minimum for the symmetric case 
Li = L2 (refer Fig. ([3])) when it equals 2~^/^+^ and tends to zero as L — )■ 00. In terms of 
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the system dimensions it is also clear that when N1/N2 is fixed and the system dimensions 
tend to infinity the skewness tends to zero. To give some numbers, for the cases shown in 
Fig. ([3]) with L = 16 and Li + L2 = 8, the average skewness 7 = .0078, .0165, .0628, and 
.2509 when Li = 4, 3, 2 and 1 respectively. These match well with the analytical estimates 
above, given that the numerical values were from 1000 trials. 



III. ENTANGLEMENT 



With the statistical properties of the partial transpose, entanglement between the sub- 
spaces 1 and 2 can now be calculated via the negativity or the logarithmic negativity. 



A. Logarithmic negativity 

The average log-negativity between two subsystems 1 and 2 of dimensions A^i and N2 
is now sought. It is assumed that the system 1+2 is a subsystem of a random pure state 
in a NiN2N-i dimensional Hilbert space, and the average is over the ensemble of uniformly 
distributed pure states in this space. Recall that the log-negativity is given by E^jy = 
log ( Ei \f^i\) = log (1-2 Ei;M,<o/^0 = log (1-2 /^<o/^ P{f^) dfi) where P{fi) is Wigner's 
semicircle given in Eq. ( fT3l) . Thus 



(Eln) = log 



- sin-i f i) + ^ W 1 - i- f 1 + 2^2 
TT \rJ 3nR\ i?2 



(31) 



where, as defined earlier, R = NR = 2^y N1N2/NS. This is valid for i? > 1. When 
R = 1 (or A^3 = 4A^iA^2), which is the critical case, this formula gives zero for the average 
log-negativity, while this is not true as discussed below. When _R < 1, PPT states are 
predominantly obtained and (Eln) = 0. Fig. @ shows how well Eq. fl3T]) works. 
On the other hand for R^ 1, deep in the NPT regime, Eq. ( 13T|) gives 



mN2 



One may compare this with the maximum possible log-negativity of a state in Hilbert space 
of dimension N1N2, as well as the average log-negativity over pure states of subsystem 
1+2. From Fig. (jl]) it is clear that there are deviations when Li + L2 = L, that is the 
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Length of block of spins (L^+L^) 



FIG. 4. The average entanglement, as measured by the log-negativity, between blocks of various 
sizes compared with the analytical results in Eq. ()3ip . The sizes Li and L2 are such that if Li + L2 
is even, they are equal and if it is odd, they differ by 1. The average log-negativity is zero when 
L1 + L2 <L/2-l. 



subsystem 1+2 is pure. If we put Ai'a = 1 (equivalently, L3 = or Li + L2 = L) we get that 
(Eln) ^ log(8v^iViiV2/(37r)). Now we present a more accurate and independent derivation 
of the average log-negativity in this case. 



1. Entanglement when pi2 is pure 

Bipartite entanglement in a random pure state is known to be very large. When ^3 = 
the Li + L2 qubits are in a pure state. The eigenvalues of are directly related to the 
eigenvalues of pi, the reduced density matrix of subsystem 1. If the eigenvalues of the latter 
are Aj, i = 1, ■ ■ ■ , 2^^, from Schmidt decomposition we have that (pi2 = I'ipu) {i^uD- 

1^12) = Yl VXl#)l#)' = E V^M^)\^?){<pf\{^?\ (33) 

i ij 

It follows that the eigenvalues of are {Aj, i^/AjAj; i 7^ j, i, j = 1, ■ ■ ■ , 2^^}, the eigen- 
vectors being and \(f)f^)\(f)^p) ±\(l)^^^)\(l)\^^) when i ^ j. The rest of the eigenvalues, 
if any, are zero. 
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Thus the average log-negativity is found as 

/ /N1N2 
(Eln) = Mog f ^ l/ij 




log 



(34) 



j=i 



Here the number k is found on using the Mar cenko- Past ur distribution of Eq. ([7]) to be 



27T .L V X 



g - 1) <^ 2F, 



1 _i2 



— 2-^1 



112 



2'2' '(VQ-l) 



(35) 



.2' 2' '(yg-i)2_ 

where x± = (1 ±l/^/Q)^, and Q = N2/N1 > 1. In the special case when Ni = N2, or Q = 1, 
the integral in Eq. ( |35|) is elementary and leads to 



(Eln) ~ log 



Ni 



(36) 



This can be compared with Eq. fl32|) which comes from a semicircle and the simple model. 
One sees that they are indeed close, but intriguingly differ by a square in the constant. We 
do not expect the semicircle to hold in the case when = 1 or Li + L2 = L. Indeed the 



density of the partial transposed spectrum is known in this case 28]. Fig. (jS]) shows the 
deviation of the spectrum from the semicircle when the system 1 + 2 is a pure state, and 
a distinctive cusp distribution is seen. The same figure also shows how poorly a semicircle 
with the same first two moments will fare. In the case when L = 6 and Li = L2 = 3 the 
scaled radius of a purported semicircle will be (from Eq. ( fT6i) ) R = 2^^/^+^^ = 16, which is 
also shown for comparison. A more detailed study of the transition from the cusp to the 
semicircle is warranted, but not carried forward here. 

We also note parenthetically that the average negativity, as defined in Eq. (jlj) when pi2 
is pure and Ni = N2 is given by 



{N{P12)) 



Ei l/^il - 1 



3tt 



(37) 



This maybe compared with an equation for the same quantity in ^] which states that 
{N{pi2)) = 0.72037(A''i — 1)/2, where the constant was arrived at numerically. One sees that 
indeed ^ = .7205..., and hence there is agreement on the principal term, while the 0{1) 
terms are however different. Indeed Eq. ( 137|) agrees very well with numerical results, the 
differences being indistinguishable from statistical fluctuations. 
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Wigner's Semicircle 
Law 




FIG. 5. Density of (scaled) /x, the eigenvalues after PT for two different cases. When L = 10 and 
Li = L2 = 3 , Wigner's semicircle law is a good fit. In the case when L = 6 and Li = L2 = 3 a 
different distribution is obtained. A corresponding semicircle of radius 16 is also shown. In general 
a semicircle is obtained for Li + L2 <^ L. 

B. Extreme Value Statistics and entanglement at the critical case 

In the critical case when = AN1N2, or in terms of the number of qubits Li + L2 = 
L/2 — 1, Eq. ( I3TI) predicts zero log- negativity as i? = 1. Thus there should be no NPT 
states. Numerical calculations however shows that there is a finite fraction of NPT states. 
Moreover, and importantly, this is not a finite-size effect. Throughout this section we assume 
the symmetric case that A^^i = N2 so that the semicircle law is valid for the eigenvalue density 
after PT. The critical case corresponds to L = 4Li + 2. See Table [T] for a calculation of 
the percentage of NPT states in several critical cases for increasing number of qubits. It 
is evident that the fraction of NPT states increases with dimensionality. While it is not 
obvious, it is argued that the fraction of NPT states saturates to a value that is close to 3% 
and 17% respectively for complex and real states. 

On the other hand when in the neighborhood of criticality (in terms of the number of 
qubits), L = 4Li + 1 and L = 4Li + 3, the percentage of NPT states rapidly increases 
to 100% and decreases to 0% respectively, see Table [Tll Also one notes that in these cases 
results have been presented for real states, as the fraction of NPT states are still significant 
for small L and the numbers are reliable. 

As mentioned earlier, for the critical case the radius of the semicircle is such that the 
hard lower limit is exactly at zero, that is the scaled radius R = 1. However it is also well 
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TABLE I. Percentage of NPT states for Li = L2 and various L for the critical case when Li +L2 
L/2 - 1. 



i 


L 


% NPT (Complex states) 


% NPT (Real states) 


1 


6 


0.06 ± 0.008 


3.18 ± 0.017 


2 


10 


1.40 ± 0.036 


7.82 ± 0.085 


3 


14 


1.92 ± 0.065 


11.18 ± 0.121 


4 


18 


2.40 ± 0.077 


13.43 ± 0.161 


5 


22 


2.60 ± 0.145 


15.17 ± 0.35 



TABLE IL Percentage of NPT states for Li = L2 and various L (Real states). 





L = 4Li + 1 


% NPT 


L = 4Li + 3 


% NPT 


1 


5 


25.39 


7 


4.4 X 10-2 


2 


9 


96.82 


11 


8.3 X 10"^ 


3 


13 


100 


15 


< 10-5 


4 


17 


^ 100 


19 


?a 


5 


21 


100 


23 


« 



known that there is a tail to the semicircle in which the extreme eigenvalues lie. The entire 
tail is then responsible for the existence of NPT states at criticality. If we are interested 
in the fraction of NPT states, this is the fraction of states such that ^mmi the minimum 
eigenvalue after FT, is less than 0. At criticality therefore it is a problem in the theory of 
extremes. In the absence of a more elaborate random matrix model, we can continue to 
use the simple model introduced earlier and see how it fares, as the theory of extremes is 
well-developed for the Gaussian ensembles. 

For N X N GUE matrices, the diagonal and off-diagonal elements ((both real and imag- 
inary parts) are drawn from the normal distributions Af{0, o"^) and Af{0, o"^ jl) respectively. 
The limits of the semicircle are ±2o"A/iV. While most of the eigenvalues lie in this range, 
some do not. The problem of estimating the number of eigenvalues outside of this range 
has been studied for long, for example see [49|]. The result about the largest eigenvalue 
distribution is now stated for cr = 1. If \max is the largest eigenvalue, then 

Xnax - 2ViV 



X 



iV-V6 



(38) 
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has a limiting distribution for large that is not one of the classical extreme value distribu 



tions, but is the Tracy- Widom distribution 
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5l| . Thus the Prob[Amaa; < a;] — -^2(2;) and 



the probability density of the scaled variable x is dF2{x)/dx, where F2{x) is obtained from 
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53l for details of a numerical 



a solution of the Painleve-II equation. See for example 
procedure that enables this. 

Applying the above to the model in Eq. (fT2|) . we need to take into account the shifted 
center and the appropriate variance of the elements of the random matrix A. Also we need 
to consider that we are interested in the minimum rather than the maximum, which is 
fixed easily as the density of the eigenvalues is symmetric about 0. Thus the density of 
the minimum eigenvalue is p{x) = —dF2{—x)/dx. Since (tr(y4^)) = N'^a'^, using Eq. ( IT5l) 
gives the variance of the diagonal elements of A to be o"^ = 1/(A^^A'3). Thus we need to 
consider (/i — 1/A^) x N^/I^ as the eigenvalue for a corresponding zero centered GUE with 
a unit variance for its diagonal elements. Thus the appropriate variable for the minimum 
eigenvalue after PT is 

X = (^^M{Nfi^,n - 1) + 2yiv) A^Ve. (39) 

As we are especially interested here in the critical case when N3 = AN = AN1N2, the variable 
X is simply y/N^N'^^^ Umin = '^N^^^ l^min- Thus the fraction of NPT states, say Jmpt is simply 
the area under the universal Tracy- Widom density, corresponding to x < keeping in mind 
that we are now dealing with the minimum eigenvalue. Thus, the simple RMT model for 
the matrix after PT results in the estimate that 

W = l-i^2(0). (40) 

Note that this is just a number (independent of matrix dimensions that are assumed to be 
large) that is numerically found to be ~ .03. For the case when L = 22, Li = L2 = 5 qubits 
we find numerically that there are 2.58% of states that are NPT, thus there is reasonable 
agreement. Figure ([6]) shows the distribution of x for two instances of critical dimensions. 
As the inset indicates, clearly there is a shift from the Tracy- Widom distribution. Indeed 
the limitations of the model of the PT as a GUE member is reflected in the statistics of the 
extremes in this way. 

One needs to add an additional shift for there to be a good match with the Tracy- 
Widom distribution. A numerically determined shift is applied to the two cases and the 
result is shown in the right panels of Figs. (|6|7|) . The shift is a positive number s such 
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The Tracy-Widom 

Distribution (Complex) 
A— A h=22, h^=h^=5 
(^L=18,L,=L,=4 




l e ' o ^ ika 



FIG. 6. Distribution of the minimum eigenvalue of various cases of Li = L2 and L. The 

right panel shows the result on using the shift which is determined numerically. This corresponds 
to complex states after rescaling as given in Eq. (|39l) . 



The Tracy-Widom 

Distribution (Real) 
A — A L=22, Lj=L,=5 
mL=18,L =r=4 




FIG. 7. Same as previous figure, but for the case of real states. 



that a; + s is given by R.H.S. of Eq. fl39|) . If the shifted distribution is used for the cases 
when (L = 18, Li = L2 = 4) and (L = 22, Li = L2 = 5), there are ^ 2.4% and ^ 2.5% of 
NPT states, in closer agreement with numerical simulations. Note that the shift will result 
in a smaller area as the right end of the integration is moved from zero to — s. This shift 
gets smaller for larger dimensionality and the fraction of NPT states seems to approach the 
fraction for the unshifted distribution. 

Whereas in this section the primary case of complex density matrices and the GUE has 
been considered, the real case is of considerable interest as well. Prevalence of time-reversal 
symmetry in many systems makes the real case important wherein the density matrix on 
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PT is modeled by matrices from the GOE. While there are no essential differences in the 
density of states after PT, both being close to the Wigner semicircle, they have very distinct 
distributions for the extreme eigenvalues. As a consequence, at the critical case, the fraction 
of NPT states is considerably higher for real states. It turns out that the correct scaling for 
the real case is the same as that in Eq. (l39l) and the distribution of the smallest eigenvalue 
after PT is shown in Fig. ([7]) for the same dimensions as for the complex case. It is clear 
from this that the fraction of NPT states is indeed larger, and is ~ 16.8%. In the real case 
too a shift is needed for good agreement with the relevant Tracy- Widom distribution, which 
is also written in terms of solutions to the Painleve-II equation 50|, |5l|- From numerical 
simulations for the cases when (L = 18, Li = L2 = 4) and [L = 22, Li = L2 = 5), the 
fraction of NPT states is 13.4% and 15.2% respectively, approaching the 16.8%. If the shift 
is incorporated and then the fraction is calculated there is good agreement even for finite L. 
Thus for example in the case when {L = 14, Li = L2 = 3) numerical simulations result in 
11.16% of NPT states, while the shift adjusted area under the Tracy- Widom density gives 
11.58%. 



1. Average log-negativity at the critical case 

Using the asymptotic Wigner semicircle results in zero log-negativity, yet there is still a 
fraction of NPT states due to the smallest eigenvalues in the tail. Although the percentage 
of NPT states can be quite high the log-negativity of the entanglement is small. Numerical 
simulations indicate that for at least the dimensions that have been considered here, if there 
are any negative eigenvalues at all on PT, there is only one. Thus the smallest eigenvalue 
/imin almost wholly controls the entanglement in the critical case. Assuming that this is the 
case gives 

ElN = log ( ^ = log (1 - 2 ^ /ii) ^ -2fXrnin 

i i^,^J.i<0 

Thus the average log- negativity at critical dimensions is given by 

2 

{Eln} ~ -2{l^rmn&{-l^min)) = ^^j^J/Q J "(^ + s)p{x)dx ~ iV~^/^ (41) 

where Q{x) is the Heaviside step function and s > is the numerically determined shift. 
The final estimate follows from the condition of criticality that = 4N1N2 = 4N . The 
Tables-Hm and [IV] show how well this estimates the average log-negativity in three cases 
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for both complex and real states, respectively. One sees that the real states have a larger 
entanglement or log-negativity in agreement with their also having a larger fraction of NPT 
states. 

TABLE III. Average log-negativity for Li = L2 and various L for the critical case (complex). 





L = ALi + 2 


Numerical (-Elat) 


{Eln) using Eq.gU 


3 


14 


7.28 X 10-6 


8.39 X 10-6 


4 


18 


9.28 X 10"'^ 


8.95 X 10"^ 


5 


22 


9.47 X 10-*^ 


9.79 X 10^*^ 



TABLE IV. Average log- negativity for Li = L2 and various L for the critical case (real). 





L = 4Li + 2 


Numerical (Eln) 


{Eln) using Eq.gl]) 


3 


14 


7.62 X 10"^ 


8.26 X 10"^ 


4 


18 


9.41 X 10~6 


9.51 X 10^6 


5 


22 


1.13 X 10-6 


1.06 X 10^6 



IV. ENTANGLEMENT AMONGST THREE COUPLED KICKED ROTORS 

To study the applicability of the results above to a dynamical system, this section studies 
a Hamiltonian system of three coupled and kicked quantum rotors or standard maps. The 
quantum standard map is one of the most important paradigmsl of quantum chaos |54| and 



has been used_extensively from early on [55| to study various phenomena such as dynamical 
localization [56]. There have been experimental realizations of the quantum standard map 
using cold atoms where dynamical localization in the momentum has been observed. Two 
coupled quantum standard maps were used to study entangling power of quantum chaos 
[s^]. More recently, there have been studies of three-dimensional (3D) kicked rotors jssi and 



many interacting kicked rotors 



59| 



A single classical standard map on the unit torus is given by the equations 

q' = q + p' (mod 1) 



K 

p = p + — sin(27rg) (mod 1) (42) 

ZTT 
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which connects phase-space variables {qi,Pi) just before two consecutive kicks which are 
separated by a unit of time. The modulo 1 conditions put the map on a phase space 
torus, which models conservative systems and Poincare surfaces of sections of two-degree- 
of-freedom systems. Much is known about the dynamics of the standard map 60l|. If K = 
then the dynamics is completely integrable. For K ^ 1 the last KAM rotational torus 
breaks, resulting in large scale diffusion in the phase space. For K < 5 the phase space is 
a mixed phase space consisting of both regular and chaotic regions. For K ^ 5, the phase 
space is nearly completely chaotic with only a possibility of finding extremely small islands 
of regular motion. 

Higher dimensional and coupled standard maps have been previously studied also because 
new phenomena such as Arnold diffusion arise 6l|. The classical coupled maps that are 
studied in this paper are given by the following canonical or symplectic transformation: 



I'i = <li+Pi (modi) 
Pi = Pi + 7r sin(27rgi 



^ sin[27r(gi + qj)] (modi) 



(43) 



j>i 



where the i^j's are parameters for respective maps and bij {j > i) are the couplings. Here 
z,j = 1,2,3 and there are three coupled rotors, with a single body potential and mutual 
couplings of two-body interactions. Higher dimensional maps such as these are only poorly 
understood. This six dimensional symplectic map is akin to Poincare surfaces of section of 
4-degree of freedom systems. However for the large parameter values that we have studied 
the maps are fully chaotic, and one may consider its quantization to be one where RMT will 
be fully applicable. The question that is being investigated is the entanglement between any 
two rotors of this tripartite system as measured by the log-negativity. 

The quantum standard map is the unitary operator corresponding to the classical map . 



It propagates states from one kick to the next. In the position representation it is 



56 



[/(n',n; K, N) 



1 



exp 



-tJM — cos 
27r 



27r 



(n + a 



exp 



ITT 

7v' 



[n 



n] 



m 

(44) 



The phase space being a torus, the quantum mechanics is on a finite dimensional Hilbert 
space of dimensionality which is related to a scaled Planck constant as N = 1/h. Thus 
the classical limit is the large A^ limit. The position kets are labeled by n = 0, A^ — 1 
with eigenvalues (n -|- a)/N. Phase-space reflection symmetry is governed by a and we 
use a = .35 below to avoid having symmetric states. The unitary operator corresponding 
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to the three coupled standard maps in Eq. (H5|) is given in the position representation by 

{n'^n'2n'^\U\nin2n^) = 

3 J- , 

H U{n[, n,- K„ Ni) J] exp | -^^V^A^^ cos 

i=l j>i ^ 

Each of the standard maps have their own dimensionahty Ni. The effective Planck constant 
is l/{NiN2N3). 

We study entanglement properties of the eigenstates of U which are the pure states of 
a tripartite system and are the stationary states as far as the quantum map is concerned. 
To be more specific, the entanglement between two rotors, 1 and 2, is studied in these 
eigenstates using log-negativity. Two parameter sets are used below 

Para. Set 1 : (K^ = 8, K2 = 7,K^ = 6, 61,2 = 1.60, 61,3 = 1.51, 62,3 = 1-42) 
Para. Set 2 : {Ki = 15, K2 = 14, = 13, 61,2 = 2.60, 61,3 = 2.51, 62,3 = 2.42). 
Such large values of K^s and 6's ensures that the individual standard maps are chaotic 
and strongly coupled with each other. Also the parameters within each set are chosen to be 
different to break any permutation symmetry effects. Using 1000 eigenstates of U the density 
of states of is shown in Fig. (IHl)) for parameter set 1, with the 64000 eigenvalues at one's 
disposal. Here we see that the distribution fits reasonably well with that of a corresponding 
Wigner's semicircle law (see Eq. ([T3]), and recall that R = 2/^/N[N^^, N = N1N2). There 
are deviations in the tail regions especially at the large eigenvalues, and agreement between 
the two distributions improves as A'"3 increases. A similar kind of behavior in the density of 
states of Pi2 for the parameter set 2 was observed and is not presented here. It is perhaps 
amusing that the Wigner semicircle appears, perhaps for the first time, in the study of a 
dynamical system, but in the properties of the eigenstates rather than the eigenvalues. 

If 4N1N2 > N3, and a typical random state is NPT, the eigenstates of the coupled 
standard maps have a log-negativity that is close to that of random states, but consistently 
slightly larger. See Fig. ([9]) for the log-negativity between maps 1 and 2 for a sample set 
of eigenstates. For parameter set 2 the average log-negativity is closer to that of random 
states compared to parameter set 1 (refer Table |V]) perhaps reflecting the increased chaos 
in the classical system, although most of the standard diagnostics of quantum chaos, such 
as the nearest neighbor spacing statistics do not differentiate between the two sets. 

The increased entanglement, as measured by the log-negativity, for the standard map in 
comparison to random states, is consistent with lowered multipartite entanglement between 



27r ' ' ^ 



N 



+ 



Ni 



(45) 
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^Nj=N,=8, N3=16 
^ N =N^=8, N3=32 
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N|J, 



FIG. 8. Density of /i, the eigenvalues after PT, from a thousand eigenstates of U, the quantum 
map of a set of three coupled standard maps. Dimensions used are A'^i = A''2 = 8 and various 
for parameter set 1, as used in the text. A vertical line at the origin has been shown to draw 
attention to the negative part of the spectrum. 

1,2, and 3, as well as with lowered entanglement between 1+2 and 3. In terms of the 
monogamy of entanglement, 1 and 2 can be more entangled as they are less entangled with 
3, as compared to a typical random state. If we view the third rotor as the environment, 
it implies a smaller decoherence for the subsystem 1+2. Thus we can say that the log- 
negativity in these cases furthers the BGS conjecture that quantum chaotic systems have 
RMT properties, but at the same time provides rather stringent and new tests for this. This 
is even more acute in the case of critical dimensions. 



TABLE V. Average log-negativity. 





N2 




Para. Set 1 


Para. Set 2 


Real random states 


8 


8 


32 


0.3567 


0.3558 


0.3491 


8 


8 


80 


0.1055 


0.1054 


0.1005 


10 


10 


10 


1.0041 


1.0035 


1.0032 


12 


12 


12 


1.0926 


1.0922 


1.0918 


14 


14 


14 


1.1678 


1.1676 


1.1669 



In the critical cases when AN1N2 = and the majority of random states are PPT, but 
there is a fraction of NPT states, the average log-negativity for the coupled standard maps 
is systematically again more than that of the random states of corresponding dimensions, 
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i 

FIG. 9. Log-negativity [Eln) for 800 eigenstates of the three coupled standard maps, the unitary 
operator li, for parameter set 1 (see text). The sohd horizontal line is the average log-negativity 
of random states of corresponding dimensions, and is given for all practical purposes by Eq. (I3ip . 



TABLE VL Average log-negativity for critical cases. 





N2 




Para. Set 1 


Para. Set 2 


Real random states 


4 


4 


64 


4.40 X 10~3 


3.90 X 10-3 


5.28 X 10""^ 


6 


6 


144 


4.85 X 10"^ 


4.35 X 10"^ 


2.24 X lO"'' 



see Table IVII The distribution of the eigenvalues after PT, near the left tail in Fig. f llUp 
differ from that of random states and are highlighted in the inset of this figure. Since the 
area under the curve for /i < is larger for the eigenstates of the coupled standard map 
than that of random states there is a larger average log-negativity for the former case. Also 
using parameter set 2 we see that the average log-negativity of the eigenstates of the coupled 
standard maps tends to that of random states as shown in Table IVI[ 

The percentage of NPT states for two critical cases is presented in Table IVIH where it 
is compared with that of real random states of corresponding dimensions. While the log- 
negativity is itself larger for the standard maps, the differences are not great. However in 
terms of the percentage of NPT states the differences between the dynamical system and the 
random states are stark. This data however does show that the percentage of NPT states of 
coupled standard map eigenstates is closer to that of real random states for parameter set 
2 than parameter set 1, again maybe a reflection of increased classical chaos, and that the 
RMT values may be reached asymptotically. However for finite quantum systems, where 
other diagnostics indicate agreement with RMT, such tests seem to show still influences of 
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a dynamical origin. 



TABLE VII. Percentage of NPT states for critical cases. 





N2 


iVs 


Para. Set 1 


Para. Set 2 


Real random states 


4 


4 


64 


29.30% 


20.37% 


7.82% 


6 


6 


144 


23.17% 


18.85% 


9.99% 



Finally the skewness of the density of states of p{l of the eigenstates of the quantum 
standard map U is compared with that of the analytical formula (Eq. fl30|) ) for the real case in 
Table f lVIII|) . One observes that as the dimension of the individual standard maps increases 
this skewness tends to that of corresponding random real states. When the dimension of 
the two standard maps 1 + 2 is small, and that of third is large the skewness approaches 
the random case, indicating once more increased decoherence from the third rotor. A more 
systematic study of the coupled standard maps, for various dynamical regimes and for 
other dimensions is postponed. The primary purpose of the present selection is to indicate 
relevant dynamical systems where we may see easily the results on entanglement of partial 
subsystems. 



TABLE VIII. Average skewness. 





N2 


iVs 


CSM Para 1 (using Eq. 


m) 


Analytical using Eq. (jSOp 


4 


4 


64 


0.1126 




6.25 X 10-2 


4 


4 


150 


8.05 X 10-2 




4.08 X 10-2 


4 


4 


200 


7 X 10-2 




3.53 X 10-2 


8 


8 


16 


0.1150 




6.25 X 10-2 


8 


8 


32 


8.14 X 10-2 




4.42 X 10-2 


8 


8 


80 


5.18 X 10-2 




2.79 X 10-2 


6 


6 


144 


5.74 X 10-2 




2.77 X 10-2 


12 


12 


12 


7.81 X 10-2 




4.81 X 10-2 


14 


14 


14 


6 X 10-2 




3.81 X 10-2 


16 


16 


16 


4.74 X 10-2 




3.12 X 10-2 


18 


18 


18 


3.86 X 10-2 




2.61 X 10-2 
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Law 




FIG. 10. Density of /x, the eigenvalues after PT, for eigenstates of the three coupled standard maps, 
for various dimensions and using the parameter set 1 (see text). Figure on the left corresponds to 
Ni = N2 = 4, = 64 and that on the right Ni = N2 = 6, N3 = 144, which are critical cases. The 
insets show an enlarged view of the region near the origin of the respective figures. A vertical line 
at the origin has been shown, as before, to draw attention to the negative part of the spectrum. 



V. SUMMARY AND CONCLUSIONS 



This paper has dealt with entanglement amongst two subsystems (say 1 and 2) of random 
tripartite pure states, using log-negativity as the measure. It is found that the state of a 
subsystem is typically NPT, and hence entangled, if the number of qubits in it (Li -|- L2) is 
larger than half the total number (L). To be precise, the number of qubits in the subsystem 
should be larger than L/2 — 1, otherwise the state is typically PPT, the critical case being 
when Li + L2 = L/2 — l. It is known that the eigenvalue distribution of the reduced density 
matrix of a subsystem is given by the Mar cenko- Past ur function; but it is found numerically 
that the same for the reduced density matrix of subsystems after PT is close to the Wigner 
semicircle law, especially when the number of qubits in the subsystems is not very close to 
the total number of qubits. A simple random matrix model, proposed herein, captures both 
the NPT-PPT transition and the spectral features after PT reasonably well. 

An analytical formula for the average log-negativity is derived using the Wigner semi- 
circle law, which is in good agreement with numerical simulations. This formula deviates 
considerably when the number of qubits in the subsystem is equal to the total number of 
qubits because the eigenvalue distribution of the reduced density matrix after PT differs from 
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the Wigner semicircle law. In this case, using tools of random matrix theory, an analytical 
expression for the average log-negativity is given that holds even if the subsystems differ in 
size. This also generalizes and augments expressions for the average negativity in 26 |. 

An exact expression for the average of the third moment of the reduced density matrix 
after PT has also been derived, both for the complex and real cases. This quantity, which is 
the first moment to deviate from that of the density matrix, is remarkable in possessing per- 
mutation symmetry amongst the three subsystems. In fact it is proved that this symmetry 
is possessed not just on average but by individual states also. Therefore it can be considered 
as a possible measure of entanglement, especially as it is also a local unitary invariant and 
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is a qudit generalization of the Kempe invariant I5 [2^ 

Using the Wigner semicircle density, the fraction of NPT states and hence the average log- 
negativity at criticality is, zero. However, a small but definite fraction of states is found to 
be NPT. This fraction and the associated entanglement is found by using the Tracy- Widom 
distribution for the extreme eigenvalues of random matrices, and is in good agreement with 
numerics, especially since it is observed that in all the cases that we have come across, 
whenever a state is NPT only one of its eigenvalues is negative. This constitutes perhaps 
the first use of the well-known Tracy- Widom distribution in quantum information theory. 

Finally eigenvectors of three coupled standard maps (or kicked rotors) were studied and 
compared with that of random real states where the parameters of the map were adjusted 
such that the classical dynamics is fully chaotic and the quantization preserves time-reversal 
symmetry. While agreeing for the most part with the results of random states, deviations are 
seen prominently at critical dimensions. The deviations are consistent with the dynamical 
states possessing marginally lower tripartite entanglement than random states. It is inter- 
esting that deviations are highlighted in quantities studied here, and hence they provide 
rather stringent tests of the BGS conjecture that random matrices are models of quantum 
chaotic systems. 

This work suggests several future directions. For example, the joint probability density 
function of eigenvalues after PT is not known, but lies presumably "between" the Laguerre 
and the Gaussian ensembles. Large deviation theory can be used to give estimates of the 
extremely small fraction of NPT states for Li + L2 < L/2 — 1 when the states are dom- 
inantly PPT. The comparison with dynamical models such as spin chains and oscillators 
will be interesting, especially in regimes where random matrix theory may not hold. Finally 
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we emphasize the occurrence of the Wigner semi-circle law in quantum chaotic systems, 
originating not in the eigenvalues, but in the properties of eigenf unctions. 



APPENDIX: AN EXACT EVALUATION OF ( tr 



Pl2 

In this appendix an exact evaluation of the ensemble average of the third moment of the 
PT, the first to depart from that of the density matrix, is calculated. Hence the skewness 
of the eigenvalue density of p\\ is found. Here the notations N-i = N = N1N2 and 
M = N1N2N3 are used. For a bipartite partition of a pure state of Li + L2 qubits and the 
remaining = L — Li — L2 qubits, a general state is given by 

N-l 

\^) = Y1 £ ain\i)\n), < Li + L2 < L. (46) 

i=0 n=0 

Hence the reduced density matrix of Li + L2 qubits {P12) and tr [(pi2)^] are given by 



{pi2)ij = 



m=0 ^^^^ 



N-l N3-I 

tr ((pi2)^) = ^ ^ aima*^ajnalnakpa^p. 



After PT (pg).. = (p^hj where 

i:= gii,j) = i-inod{i,N2) + mod{j,N2), 

(48) 

j ■= gU, = J - mod(j, A^2) + mod(z, N2). 
The function = {g{h j), g{j, i)) is bijective and is its own inverse, since perform- 

ing PT twice keeps elements of pi2 unchanged. This implies that i = g{i,j) and j = g{j,i), 
i.e. (pfl).. = ipuYfi- Thus we state the following simple conclusions are useful lemmas: 

Lemma li = jiSi = j. 

Lemma 2 g{i,i) = i. 

Lemma 3 If ^ (gii, j), gij,i)) and ^ (gii' J'), gif and j 7^ j' {i 7^ i'), 

then g{i,j) 7^ g{i'-,j') {g{jii) 7^ gii'i^'))- In words, elements that differ in column (row) 
position get mapped after PT to positions that differ in row (column). 

The expression for tr[(p^2)'^] using Eqs. (H7|) and the above function is 

N~l Nz-l 
i,j,k=0 m,n,p=0 



32 



where h = g{i,j), ji = gij,i), 32 = 9{j,k), h = g{k,j), h = g{k,i) and is = g{i,k). 
Using this one obtains the ensemble average (tr[(p^2)^]) follows. Here the fact that after 
averaging only even powered terms will be nonzero and odd powered terms will be zero is 
used. It can be seen in Eq. ( H9|) that there are three possible cases for m, n and p, namely 
m^n^p,m = n^p (which is the same asm ^ n = p and m ^ p = n) and m = n = p. In 
each of these cases the number of terms that do not vanish on averaging is first calculated, 
the last case requiring a somewhat detailed analysis. 

Case: m 7^ n 7^ p 

In this case, the only possible non-vanishing term after averaging is a product of three 
unequal quadratic terms i.e. ii = ji, j2 = ^2 and k^ = i^, which implies that g{i,j) = g{j, i), 
9{jjk) = g{k,j) and g{k,i) = g{i,k). As a consequence of the Lemma 1 above, one then 
obtains that i = j = k. Thus in this case there are exactly N non-vanishing terms on 
averaging Eq. fH9|) . each of the form |ajmp|aj„p|ajpp. 

Case: m = n 7^ p 

Now there are two possible types of non- vanishing terms, one a product of three unequal 
quadratic terms and the other a product of one quadratic and one quartic term. In the 
former possibility, h = h, k = h, ji = 32 and ii ^ ji i.e. g{k,i) = g{i, k), g{i,j) = g{k,j), 
9(3^ = dU^ k) and g{i,j) 7^ g{j, i). Again using the Lemma 1, this implies that i = k ^ 3. 
Thus there are A^(A^— 1) non-vanishing terms on averaging, each of the form |aimPl'2jmP|ctipP- 

In the later case, k^ = i^ and ii = ji = 32 = k2 i.e. g{k, i) = g{i, k) and g{i,j) = g{j, i) = 
g{j, k) = g{k,j) which implies i = k = j. Thus there are non- vanishing terms of the form 
la- l^la- P 

Case: m = n = p 

In this case there are three possible types of non-vanishing terms: (1) one sextic term, (2) 
one quadratic and one quartic term and (3) three unequal quadratic terms. 

1. One sextic term. This occurs when ii = ji = 32 = k2 = k^ = is, that is g{i,j) = 
g{j, i) = g{j, k) = g{k,j) = g{k, i) = g{i, k). This in turn implies that i = j = k. Thus 
there are A^ non- vanishing terms of the form |aim|^- 
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2. One quadratic and one quartic term. Corresponding to such terms there are two cases. 

(a) ii = k2, ji = j2 = h = i3 and ii 7^ ji. 

In terms of function g this condition is g{i,j) = g{k,j), g{j,i) = g{j,k) = 
g{k,i) = g{i, k) and g{i,j) 7^ g{j,i)- This imphes, using the Lemmas above, that 
i = k and g{j,k) = k but k 7^ j. Thus j — mod{j,N2) + mod(fc, A^2) = k, i.e. 
j — mod(j, A^2) = k — mod(A;, A^2) and k 7^ j. In other words, j and k have to be 
distinct, but have the same quotient on division by A^2- Here j and k takes values 
from to — 1 where = NiN2. To find the number of {j, k) pairs that satisfy 
these conditions, imagine dividing into A^^i intervals each of length A'"2. If j 
is selected from the possible numbers, this also fixes one such interval. The 
number k must necessarily be in this interval, but must not be j, which gives a 
choice of multiplicity (A'2 — l).Thus there are N[N2 — 1) non- vanishing terms of 
the form |ajm|'*|cffcmP- 

(b) ii = is, ji = j2 = k2 = h and ii ^ ji. 

In terms of the function g this condition is g{i,j) = g{i,k), g{j,i) = g{j,k) = 
g{k,j) = g{k,i) and g{i,j) 7^ g{jii)- Using the Lemmas above implies that 
j = k and g{j,i) = j but i 7^ j. Thus j — mod(j, A^2) + mod(2, A^2) = j i-e. 
mod(j, A^2) = mod(2,A^2) and i 7^ j. In words, one must count the number of 
distinct pairs of integers having the same remainder on division by A^2- This is 
easily seen from an argument similar to that in the previous paragraph to be 
N{Ni — 1), each corresponding to a non-vanishing term of the form |ajmP|ojrra|*^- 

Combining these two cases and including the cyclic permutation of k), the total 
number of non- vanishing terms of this kind are 3A^(A'^i + A^2 — 2). 

We note in parenthesis that this is the first instance that the counting is different from 
that for evaluating (tr[(pi2)'^]) in which one has the condition i = j ^ k and its cyclic 
permutations, implying 3N{N — 1) non-vanishing terms. 

3. Three unequal quadratic terms. For three unequal quadratic terms there are fifteen 
possible cases, of which only seven are distinct up to cyclic permutation of indices 
{i,j,k). These are listed in Table (jlXp . and subsequently analyzed. 

(a) {ii = ji) ^ {j2 = h) ^ {h = h)- 
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TABLE IX. 



(a) 




= 3\) + ih 




= h) 


(b) 


{h 


= h) + ih 


= h) + 


= ^3) 


(c) 


{h 


= ji) ¥= {h 


= ■is) 7^ 


= h) 


(d) 


in 


= 32)^ (ii 


= h) 7^ (^2 


= k) 


(e) 


in 


= h)i^ ih 


= ^3) 7^ (^2 


= h) 


(f) 


in 


= h) + (ii 


= h) + ih 


= n) 


(g) 


ih 


= h) + (ji 


= h) + i^2 


= h) 



In terms of function g this condition is igii,j) = gij,i)) 7^ igij,k) = gik,j)) 7^ 
igik, i) — gii, k)). This gives i — j,j — k, i — k, j, j ^ k and i ^ k, which 
are incompatible conditions and hence the number of terms of this kind is zero. 
The multiphcity of this case under cychc permutation of the labels k) is 1. 

ih = ji) 7^ (j2 = h) 7^ ih = ^3)• 

In terms of function g this condition is igii,j) — gij,i)) 7^ idiJ^k) — gik,i)) ^ 
(fl'(^)j) — 9ihk))- The first equality gives i — j, and the inequation gij,k) 7^ 
gik,j) implies that j 7^ k. However the second equality (using i — j) implies that 
i = k. Thus these conditions are incompatible and the number of terms of this 
kind is zero. The multiplicity of this case under cyclic permutation of the labels 
iij,k) is 3. 

ih = ji) 7^ O2 = ^a) 7^ (^2 = ^~3)- 

In terms of function g this condition is igii,j) = gij,i)) 7^ igij,k) = gii,k)) 7^ 
idi^yj) = 9ik^i)). This first equality gives i = j, while the inequation gii,k) 7^ 
gik,i) implies that i ^ k. Further gii,j) = gii,i) — i ^ oih^) ^ mod(i,A^2) 7^ 
mod(/c, A^2) and i 7^ gik, i) =^ i — mod(i, A^2) k — mod(/c, A^2) i-e. remainder 
and quotient of i and k, under division by A^2, are not same. Again dividing N 
into A^i intervals each of length N2 gives NiN — Ni — N2 + I) as the number of 
possible triples (z, j, k) that satisfy these conditions. The multiplicity of this case 
under cyclic permutation of the labels (i, j, k) is 3. 

ih = J2) 7^ (ji = ^3) 7^ ih = is)- 

In terms of function g this condition is igii,j) = gij,k)) 7^ igijji) = gik-ii)) 7^ 
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(fl'(^)j) — 9{hk))- The inequations and the Lemma 1 above imply that i ^ 
j ^ k. It is then not hard to see that the equahties are incompatible with this 
condition. For instance (gii, j), 9ij,i)) whereas {j,k) {g{j, k), g{k, j)), 

and using Lemma 3 with = {j,k), g{i,j) 7^ g{j,k), which violates one of 

the requirements. Thus the number of terms of this kind is zero. The multiplicity 
of this case under cyclic permutation of the labels k) is 3. 

(e) {k = 32) 7^ (ii = ^3) 7^ (fc~2 = ^~3)- 

In terms of function g this condition is {g{i,j) = g{j,k)) ^ {g^j-ii) = g{hk)) 7^ 
{gik-ij) — g{k, i))- An analysis very similar to the previous case shows that these 
conditions are incompatible too. The multiplicity of this case under cyclic per- 
mutation of the labels k) is 3. 

(f) {H = k2)j^{jl=h)^{j2=(3). 

In terms of function g this condition is {g{i,j) — g{k,j)) ^ igij,i) — g{k,i)) 7^ 
{g{j,k) — g{i,k)). The inequations again imply that i j k, however the 
Lemma's do not lead to incompatible conditions. 

g{h j) = g{k, j) ^i- mod(i, N2) ^ k - mod(/c, N2) 
gU, i) = g{k, i)^ 3- mod(j, N2) ^ k - mod{k, N2) 
g{3: k) = g{i, k)^i- mod(i, A^2) = j - mod(j, N2). 

Thus i, j, k have the same quotient on division by N2, and as they are all distinct 
they have different remainders. To find the number of triples k) that satisfy 
these conditions, once more divide an interval of length into A^i intervals of 
length A^2- One can select i in N possible ways, which fixes the quotient on 
division by N2- The integers 3 and k must then be one of the possible N2 — 1 
numbers, without also being equal. Hence the number of terms of this kind is 
A'"(A^2 — 1)(A^2 — 2). The multiplicity of this case under cyclic permutation of the 
labels /c) is 1. 

(g) («i = ^3) 7^ (ii = 32) 7^ {k2 = h). 

In terms of function g this condition is {g{i,j) — g{i,k)) ^ {g{j,i) — gij,k)) 7^ 
{g{k, j) — g{k,i)). A similar analysis as for the previous case shows that i, j 
and k are distinct but they have a common remainder on division by Thus 
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there will be N{Ni — l)(A^i — 2) number of non- vanishing terms in this case. The 
multiplicity of this case under cyclic permutation of the labels k) is 1. 

Thus the number of terms with a product of three distinct quadratics, including the 
multiplicities is ?,N{N -Ni-N2 + l) + N{Ni - l){Ni - 2) + N{N2 - 1){N2 - 2). In 
contrast the number of such terms in the evaluation of (tr[(pi2)^]) is N{N — 1){N — 2). 
In general note that one can recover results for the density matrix prior to PT from 
those after PT by replacing N2 ^ 1 and A^i — )■ A^. Thus the results for PT present a 
particular generalization. 



631 ] for the case of complex states are stated 



The exact RMT ensemble average values 
below for convenience: 

6 2 

^'""'"''^ " M(M+l)(M + 2)' = {\a^n\'\aim\') = M{M + l){M + 2y 

(lo-inl I'^jml ) (lo-ml |Q-jm,| |Q-A;p| ) (lo-ml I'^jm] \(^ip\ ) (I'^inl |'^jn| l^-fepl ) 

n |2| |2| |2\ _ n |2| |2| |2\ _ /i |2| |2| |2\ _ \ 

VPinl |%n| IQ-jpl / — \|Om| Pjra| Pjpl / — \\Clin\ |Q-jn| |Q-A;ra| / — _j_ \"^(^]\/[ _|_ 2) ' 

These averages are multiplied by the number of non-vanishing terms and by the respective 
multiplicity for m, n and p (for m 7^ n 7^ p it is N^{N^ — 1){N^ — 2), for m = n ^ p ii 
is 3A^3(A'^3 — 1) and for m = n = p it is simply N^), and added together. Use is made of 
A^A^3 = N1N2N3 = M and a straightforward simplification of 22 terms results in significant 
cancellations, leaving just 4 terms finally. This results in 

.^ r^.s. _ Nl + iV| + Nj + 3N,N2N, 

with the remarkable permutation symmetry evidently displayed. 

A similar analysis can be done for the case of averaging only over real states. The counting 
remains identical to the complex case, while the averages differ as: 

/ 6 \ /42\ /42\ 

^"'"^ " M(M + 2)(M + 4)' ^"^^^^J"/ = = M(M + 2)(M + 4)' 

/ 2 2 2 \ / 2 2 2 \ / 2 2 2 \ / 2 2 2 \ 

X^in^jm^kp) ~ X'^in^im^kp) ~ X^in^im^ip) ~ X^in^jn^kp) ~ 

/222\_/222\_/222\ ^ 
X^in^jn^ip) ~ X^in^jn^jp) ~ X^in^jn^kn) 



M(M + 2)(M + 4) 

This leads to the ensemble average: 

/. ^ T,^3, _ Nf + Nj + Nj + 3{N, + N2 + N, + N,N2Ns) 

^^^^^^ ^ " {N,N2N, + 2){N,N2N, + 4) ' ^^'^ 
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where we have used M = N1N2N3. 

Using the earher statement that the averages prior to PT can be found from those after 

by the replacement A^2 — ^ 1 and Ni ^ N = N1N2, one gets the 

, . .3 ^ NfNl + iV| + 3 + 1 

^ '^^''^ ^ {N,N2Ns + l){mN2N, + 2) ' 

^ '^^''^ ' (iViiVsiVs + 2){N,N2N, + 4) 

for the case of complex and real cases respectively. This indeed agrees with a previous 
calculation of this quantity, Eq. (5.11) in 32|, where the complex case is considered. 
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